Scattering amplitudes at finite temperature 
M.E. Carrington"'^ Hou Defu'^-'' and R. Kobes^-'f] 

" Department of Physics, Brandon University, Brandon, Manitoba, R7A 6A9 Canada 
^ Winnipeg Institute for Theoretical Physics, Winnipeg, Manitoba 
Institute of Particle Physics, Huazhong Normal University, 430070 Wuhan, China 
University of Winnipeg, Winnipeg, Manitoba, R3B 2E9 Canada 

We present a simple set of rules for obtaining the imaginary part of a self energy diagram at finite 
temperature in terms of diagrams that correspond to physical scattering amplitudes. 

PACS numbers: PACS: ll.lO.Wx, 11.15.TK, 11.55.Fv 



I. INTRODUCTION 

In this paper we discuss a set of rules for calculating the imaginary part of self energy diagrams. These rules have a 
simple diagrammatic representation in terms of scattering amplitudes. They have been deduced by studying a large 
set of diagrams; a derivation from first principles is in progress and will be presented in another paper. 

It is well known that the imaginary part of the retarded self energy is an important quantity in thermal field 
theory: it provides information about decay and production rates of particles, among other things. The physics 
that is contained in the imaginary part of the self energy is revealed by writing it as the product of two scattering 
amplitudes. At the one loop level, the structure of the scattering amplitude has been understood for some time [Q. 
Extension to higher loops is not straightforward. In this paper we discuss cutting rules in the context of this problem. 
We show that rules exist that make it easy to understand the physical content of the imaginary part of a self energy 
diagram at high loop order. 

We begin by reviewing some basic concepts and defining some notation: 

a) Off shell propagators describe the behaviour of fields that propagate through the medium until they undergo 
interactions with other fields. Diagrammatically, off shell propagators are represented by continuous lines that begin 
and end at interaction vertices. 

b) On shell propagators carry delta functions of the form 5{P^ — m?) where P is the momentum of the field. On shell 
propagators correspond to fields that do not propagate through the medium. 

c) Consider any closed loop in which all of the propagators are off shell and momentum is free to flow around the 
loop. If any one of the propagators in the loop is put on shell, momentum is no longer free to flow around the loop 
and the loop is effectively 'opened.' 

d) There are two kinds of propagators that are on shell. We call these two types of propagators "cut" propagators 
and "tic-ed" propagators. Cut propagators and tic-ed propagators carry different thermal factors. This point will be 
explained in detail. 

e) A "cut line" is a line that divides the self energy into two pieces, each of which has one external leg. Any propagator 
that is crossed by a cut line is put on shell and becomes a "cut" propagator. 

f) Diagrammatically, our notation is as follows. In a self energy diagram, a cut propagator is a propagator that is 
crossed by the cut line and a tic-ed propagator is drawn with a double tic mark through it. To obtain scattering 
amplitudes, all on shell propagators (cut or tic-ed) are split into two pieces, each of which has an end that is not 
connected to a vertex. As a result, the cut line divides the self energy diagram into two separate amplitudes, and 
(as will be explained below) the tic-ed propagators cause each amplitude to have the form of a tree amplitude, with 
no closed loops. The lines obtained from the splitting of on shell propagators represent the emission or absorption 
of fields by the medium. The lines that represent absorbed fields are drawn slanting backwards from the vertex and 
lines that represent emitted particles slanting forward from the vertex. 

g) For any given diagram, the number of 5 functions, or the number of on shell = (cut + tic-ed) propagators, is equal 
to L + 1 where L is the number of loops. 

To begin, we consider a one loop calculation. Following Weldon we look at the simple case of a scalar field $ 
coupled to two other scalars 4)i and 4)2 through a cubic interaction. The production rate for the field $ is obtained 
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from the imaginary part of the one loop self energy shown in Fig. [la]. We obtain the imaginary part by drawing a 
cut line through the diagram. At one loop, there is only one way to draw a cut line so that each half of the diagram 
contains one of the two external legs (Fig. [lb]). This cut line produces two cut propagators, each of which carries a 
delta fmiction (there are no tic-ed propagators in this case). 




(a) (b) 



FIG. 1; One- loop self-energy 

Expanding these two delta functions we find that the imaginary part of this diagram can be written as the sum of 
four terms, each of which corresponds to the square of a scattering amplitude j^]. We obtain, 

2ti 

— — [5{qo - t^i - t^2)[(l + "i)(l + 7^2) - nin2\ 
'1ZUJ2 

+ S{qo + uji - LL!2)[{ni{l + 712) - (1 + ni)n2] 
+ 5{qo - uji + t^2)[(l + ni)n2 - ni(l + 712)] 

+ 5{qo +UJI+ UJ2)[ni722 - (1 + "i)(l + "2)]] (1) 

where oji = yj^ + mf , uj2 = \/{p + O)^ + ^2 ni = 1/ (e''"^ — 1) with (3 the inverse temperature and g the coupling 
constant. The first term in this expression corresponds to the probability for the decay <& — > 0102 with a statistical 
weight (l-|-ni)(l-|-n2) for stimulated emission, minus the probability for the inverse decay 0i02 with the weight 
nin2 for absorption. Note that the thermal factor (1 -I- ni)(l -t- 712) — ^1^2 which reflects the physics of the process 
involved could be written in the mathematically simpler form 1 -(- ni -I- ^2 at the cost of loosing information about 
the physics. Similarly, the second term gives the probability for the decay <i>(/)i — > 02 (which involves the absorption 
of a 01 field and the emission of a 02 field), minus the probability for the reverse process 02 — >■ <l'0i, with appropriate 
thermal weights. The interpretation of the third and fourth terms is straightforward. These four process are shown 
in Fig. [2]. 



/d' 
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FIG. 2: The scattering amplitudes corresponding to Fig.||l|] 

We would like to study what happens at higher loops where scattering amplitudes have a much more complicated 
structure. The task of calculating the imaginary part of the self energy and separating it into scattering amplitudes 
is complicated in different ways, depending on the technique that is used to calculate the imaginary part of the self 
energy. 

In the imaginary time formalism (ITF), one starts in Euclidean space, calculates Green functions with imaginary 
time arguments, and performs an analytic continuation to real time at the end of the calculation. One attractive 
feature of the ITF is that it satisfies the intuitive belief that it should be possible to obtain finite temperature field 
theory from zero temperature field theory by adding thermal weights to the Feynman rules, in some fashion. The 
ITF is mathematically simpler than the RTF because of this structure. However, the price one pays for mathematical 
simplicity is that some physical processes are hidden. Thinking more carefully we realize that this feature of the 
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ITF is not surprising since we should not in fact expect that finite temperature field theory should have the same 
structure as zero temperature field theory. In a thermal situation, individual fields do not simply propagate through 
the vacuum, but interact with a medium. Consequently, a specific scattering amplitude will involve a number of 
interference process that are not present at zero temperature. All of these processes are, of course, present in the ITF 
calculation, but the compactness of the notation effectively hides them from view. A lengthy procedure for extracting 
physical amplitudes from the ITF has been discussed by Wong j^. This extremely complicated calculation has been 
carried out explicitly by Majumder and Gale for the two loop vector boson self energy in thermal QCD 

In the real time formalism (RTF), one works in Minkowski space and Green functions with real time arguments are 
obtained directly. It is well known that the RTF is mathematically more complicated than the ITF because of the 
doubling of field degrees of freedom. In the closed time path (CTP) representation of the RTF, the contour has two 
branches: the top one (Ci) runs from negative infinity to positive infinity, and the bottom one (C2) runs backwards 
in the other direction. These two branches give the propagator a matrix structure. The four components of the 2x2 
matrix are labeled -D12, £'21, and D22 and correspond to propagation along Ci, propagation from Ci to C2, etc. 
The propagator Dn corresponds to time ordered propagation. The advantage of working in the RTF is that is it 
easier to separate the imaginary part of a self energy into physical scattering amplitudes. This point will be discussed 
in detail. 

A great deal of work has already been done on the development of rules for the calculation of the imaginary parts 
of diagrams in the RTF. One set of rules was derived by Kobes and Semenoff [Q using the 1/2 representation of the 
RTF. A set of rules based on more general representations, such as the R/A or Keldysh representaions, has been 
developed by Gelis [||. We have checked that these rules and the ones discussed in this paper are equivalent, as they 
must be. However, some intricate cancellations are required to extract the scattering amplitudes from the rules of 
[Q, ||]. Several other authors have looked directly at scattering amplitudes. Aurenche and collaborators have separated 
the scattering amplitudes contained the imaginary parts of the two loop photon self energy diagrams by dividing the 
phase space of the momentum integrals into regions that contain the different possible combinations of signs of the 
frequencies of the fields However, this technique would be difficult to generalize at higher loops. Brandt and 
collaborators have developed a diagrammatic representation for retarded green functions in terms of tree scattering 
amplitudes in the high temperature limit Q. 

The paper is organized as follows. In section II we present the first part of our cutting rules: how to determine 
the diagrams that contribute to the imaginary part of a given self energy, and for each diagram, which propagators 
are on shell, retarded, or advanced. Section III contains an discussion of how to interpret the diagrams produced by 
our rules in terms of scattering amplitudes. Section IV contains the second part of our rules: how to determine the 
thermal factor for each propagator. Section V contains a list of the self energy diagrams that were used to deduce 
these rules, and the diagrams that contribute to their imaginary parts. In section VI we present our conclusions. 



II. FIRST PART OF THE CUTTING RULES: PROPAGATORS 



A. Allowed Diagrams 



For any self energy diagram, draw all possible cut lines (a cut line is any line that divides the diagram into two 
pieces, each of which contains an external leg). A cut line that opens all loops will produce L + 1 cut propagators and 
thus L + 1 delta functions. No other propagators will be on shell. Other cut lines will leave some loops unopened, and 
produce less than the required number of L + 1 delta functions. Add tic marks to these diagrams in every possible 
way so that 1) all loops are open; 2) the total number of delta functions is equal to the L + 1 and 3) it is possible to 
move from either external leg, to the cut line by following a continuous path along uncut and untic-ed propagators. 
As an example, for the scalar self energy in Fig. [3], the allowed cut diagrams are shown in Fig. [4] (for a scalar theory, 
the diagrams in Figs. [4a, b] and Figs. [4c, d] are equivalent). 



FIG. 3: A two-loop self energy 




(d) (e) 



FIG. 4: The cuts that contribute to the imaginary part of Fig. [3] 

Two examples of diagrams that violate rule 3) are shown in Fig. [5]. Both of these diagrams contain the right number 

of delta functions (5), and every loop is opened, but it isn't possible to get from the right leg of cither diagram to the 
cut line without going through a cut or tic-ed propagator. Neither of these diagrams should be drawn. 




FIG. 5: Two diagrams that violate rule 3) 



B. Uncut and Untic-ed Propagators 

Propagators that are uncut and untic-ed are either retarded or advanced. There is a rule to determine which. 

Consider the imaginary part of a retarded self energy. Start from cither external leg and trace any continuous path 
through the diagram (continuous means not traveling along a propagator that is cut or tic-ed) . If the momentum of a 
given propagator flows in the same direction as the path you are tracing, the propagator is retarded. If the momentum 
flows in the direction opposite to the direction of the path, the propagator is advanced. This rule is illustrated in 
Fig. [6]. 
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FIG. 6: An illustration of the propagator rule. 



III. SCATTERING AMPLITUDES 



In this section, we discuss how to interpret the diagrams produced by our rules as scattering amplitudes. 



A. QED at One Loop 



Throughout most of this paper, we will work with theory. This toy theory allows us to avoid a lot of the 
mathematical complexities associated with physical theories which are not directly relevant to the problem of how 
to calculate the imaginary part of a self energy, and how to interpret the result in terms of scattering amplitudes. 
However, in this section and the next, in order to discuss the physics of scattering amplitudes in more familiar terms, 
we will switch to QED. The structure of the scattering amplitudes in QED and (j)^ theory is similar since both theories 
have an interaction involving three fields: (j)^ theory has a cubic interaction, and QED has a photon-electron-positron 
interaction. In the diagrams in this section and the next, all arrows give the direction of flow of lepton number. We 
will look at the photon self energy. 

We start with the one loop diagram (Fig. [7a]). There is only one possible cut line that separates the diagram into 
two pieces, so that each contains an external leg. The amplitudes obtained from the right hand side are shown in 
Fig. [7b] and have the same structure as those shown in Fig. [2] . As discussed in the introduction, a fermion line that 
represents an emitted (absorbed) fermion is c;onnected to a vertex at only one end and is drawn slanting forwards 
(backwards) from that vertex. Fig. [7b] shows the amplitude for a photon to decay into an electron-positron pair 
(7 e+e"), the amplitude for a photon to absorb a positron from the medium and emit an electron (76+ — > e~), 
etc. In the future, we restrict to positive frequencies for the external field. In this case the process represented by the 
last amplitude is kinematically forbidden. The left hand side of the cut self energy in Fig. [7a] gives the conjugate of 
the amplitudes in Fig. [7b] , and thus the product of the right hand side and the left hand side is a real probability. 




(a) 




FIG. 7: One loop photon self energy and scattering amplitudes. 
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B. QED at Two Loops 

There are two diagrams that contribute to the photon self energy at two loops. They are shown in Figs. [8a, b]. For 
each of these diagrams, there is more than one possible cut line. 




(a) (b) 



FIG. 8: Two loop photon self energies. 



1. Central Cuts 

To begin with, for both diagrams, we look at the cut line that goes directly through the center as shown in Fig. [9]. 




(a) (b) 



FIG. 9: Central cuts. 

In both cases, all loops are opened and thus there are no tic-ed propagators. Both cut lines cross one photon 
propagator and two fermion propagators. The forward scattering amplitude is given by the terms in which all three 
particles are emitted. The Compton scattering amplitude occurs when the photon is emitted, and one fermion is 
emitted and one is absorbed. The amplitude for pair production is produced when the photon is absorbed and 
both fermions are emitted. When we restrict to positive external frequencies, no other possibilities are kinematically 
allowed. We describe these processes below. 

1) Forward Scattering 

The amplitude produced by left hand side of Fig. [9a] is shown in Fig. [10a]. The right hand side of the diagram is 
the conjugate amplitude. There is another contribution to the self energy that is the same as Fig. [8a] except with the 
propagator correction on the lower line, or alternatively, with the flow of lepton number in the fermion loop routed in 
the opposite direction. This graph is just the hermitian transpose of Fig. [8a] and the central cut produces the same 
amplitude with the roles of the electron and positron reversed, as shown in Fig. [10b]. Again, the right hand side of 
the diagram is the conjugate amplitude. These two diagrams represent the amplitude for a photon to decay into an 
electron-positron pair where one of the fermions emits an additional photon. 




The amplitude corresponding to the left hand side of Fig. [9b] is also shown in Fig. [10a]. In this case however, the 
right hand side of the diagram does not give the conjugate amplitude. Instead, we get the cross product of Fig. [10a] 
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with the conjugate of Fig. [10b]. By taking the cut Hne diagonally in the opposite direction, we obtain the reversed 
cross product, with the electron and positron switched. Combining these results we obtain the square of the amplitude 
shown in Fig. [11], which is proportional to the photon decay probability. 




FIG. 11: Photon decay. 

2) Compton Scattering 

We look at exactly the same cuts as above, but consider the amplitudes that correspond to one emitted photon, 
one emitted fermion and one absorbed fermion. The graph in Fig. [9a] produces the squares of the two amplitudes 
shown in Fig. [12a]. The diagram with the propagator correction on the bottom line gives the same amplitudes, with 
the direction of Icpton flow reversed. The diagram in Fig. [9b] gives the two cross terms shown in Fig. [12b] . 




FIG. 12: Compton scattering amplitudes from Fig. \ \L1\ . 

The graphs that are the same as those in Fig. [12b] but with the fermion lines reversed are produced by the cut 
on the opposite diagonal. Combining these results we obtain the square of the amplitude shown in Fig. [13], which is 
proportional to the Compton scattering rate. 




FIG. 13: Compton scattering. 
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3) Pair Production 

Following the same method, we obtain the square of the amplitudes for pair production as shown in Fig. [14]. As 
before, the cross term is produced by Fig. [9b] . 







FIG. 14: Pair production. 



2. Non-Central Cuts 



From this point on wc will discuss only forward scattering amplitudes. Also, we will not discuss the fact that each 
amplitude has an partner in which the flow of fermion number is reversed. These graphs are obtained in exactly the 
same way as was described above for the central cuts. 

For both of the diagrams in Fig. [8] there is a cut that does not open both loops. For the diagram in Fig. [8a] this 
cut leaves the loop formed by the propagator correction unopened. There are two ways to open this loop with a tic. 
These two graphs are shown in Fig. [15]. 





FIG. 15: Non-central cuts for Fig. [8a]. 



The amplitudes corresponding to the left hand side of Fig. [15a] are shown in Fig.[16a,b]. The first of these graphs 
represents the amplitude for a photon to decay into an electron-positron pair, where the electron subsequently absorbs 
a photon from the medium and emits it back into the medium. The second graph is the same as the first, except 
the photon is emitted into the medium and then absorbed. This photon is sometimes referred to as a spectator field. 
Spectators are on shell particles from the heat bath that enter with the in-state and leave with the out-state without 
having interacted with the rest of the participants. The right hand side of Fig. [15a] is shown in Fig. [16c]. The lone 
wiggly line represents the spectator particle from the heat bath. This spectator field is normally not drawn, but its 
presence is necessary to obtain an interference effect: when calculating the convolution of the two amplitudes, it is 
necessary to have the same number of incoming and outgoing particles. We can write a schematic equation to describe 
these amplitudes: 

(77 e"'"e~7) (g) (7e"'"e~ ; 7^7) 



where the tilde indicates the spectator field. 



(a) (b) (c) 

FIG. 16: Amplitudes from Fig. ||a]. 

The amplitudes from the left hand side of Fig. [15b] are shown in Fig.[17a-c]. Fig. [17a] represents the amplitude for 
a photon to decay to an electron-positron pair, where the electron subsequently absorbs a positron from the medium 
and emits it back into the medium. In Fig. [17b] this spectator positron is emitted and then absorbed. Fig. [17c] shows 
the amplitude from the right hand side of the self energy, with the spectator positron. Schematically we write, 

(76^ e+e^e^) (7 — > e+e~ ; — > e'*') 




(a) (b) (c) 



FIG. 17; Amplitudes from Fig. [^]. 

A real probability is obtained by combining with the amplitudes produced by the graphs that are just like those 
in Fig. [15], but with the cut line on the left hand side. All of these graphs have been extracted from the imaginary 
time calculation by the authors of ||] and In [Q they appear as Figs. [16, 17] and in ||] they are Fig. [12]. In the 
imaginary time calculation, a great deal of tedious effort is required to separate these physical contributions. Using 
our rules however, they appear immediately. 

Now we consider the non-central cut for Fig. [8b]. The unopened loop can be opened with a tic on the photon line 
Fig. [18a], or a tic on the fermion line Fig. [18b] (tic-ing the upper fermion line would change the spectator fermion to 
a spectator anti- fermion) . The amplitudes that result from Fig. [18a] are shown in Fig.[19a-c]. The first two graphs 
show the decay of a photon into a electron-positron pair with a spectator photon. In Fig. [19a] the electron absorbs 
a photon from the medium and the positron emits a photon, and in Fig. [19b] the electron emits a photon and the 
positron absorbs a photon. Fig. [19c] shows the amplitude from the right hand side of the self energy in Fig. [18a], with 
the spectator photon. We write schematically, 

(77 — > e^e~7) (g) (7 ^ e~^e~ ; 7^7) 

The amplitudes from Fig. [18b] are shown in Figs. [20a-c] . Schematically we write, 

(76" e^e~e~) (7 ^ e'^e" ; ^ e~) 




(a) (b) 



FIG. 18: Non-central cuts for Fig. [8b]. 
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(a) (b) (c) 

FIG. 19: Amplitudes from Fig. ||a]. 




(c) 

FIG. 20; Amplitudes from Fig. [^]. 
These diagrams are Figs. [19,20] in ^ and Figs. [10,11] in ||^. 

C. (j)^ Theory at More Than Two Loops 

In order to discuss scattering amplitudes at higher loop levels, we will simplify the notation as follows. 

1) We revert to (f)^ theory. 

2) We will only draw the forward scattering amplitude (all cut fields on the positive mass shell) associated with the 
left hand side of a given cut self energy. 

3) We simplify the diagrammatic representation of tic-ed fields as follows. All tic-ed propagators correspond to 
fields that are emitted and then absorbed, or absorbed and then emitted. Thus, if there are two tic-ed propagators 
on the left hand side of a cut self energy, there are four different scattering amplitudes corresponding to the four 
different ways in which these two fields could be emitted and absorbed from the heat bath. Using the notation from 
the first part of this section we would draw four diagrams, one for each of these processes, with the tic-ed propagators 
represented by lines slanting forward or backwards to represent emission or absorption respectively. In this subsection, 
we will draw only one diagram instead of four, and represent emitted or absorbed fields by arrows pointing straight 
up, or slightly slanted to one side or the other in the case where two fields are emitted/absorbed from the same vertex 
and could not be seen on the figure if they were drawn on top of each other. 

4) All propagators will be drawn as lines slanting up to the right. Note that if we define momentum variables 
so that momentum fiows to the right, all propagators are retarded (following the rule in section IIB), and thus the 
scattering amplitude has a causal interpretation: time flows to the right along any continuous line. 

5) Cut propagators represent fields that are emitted and absorbed on opposite sides of the diagram. These cut lines 
are drawn as horizontal arrows pointing to the right. 

An example of this notation is shown in Fig. [21]. In the next section we describe in detail how to determine the 
appropriate thermal factor for a given cut diagram by looking at the corresponding scattering amplitude. 
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FIG. 21: An example of the notation to be used in section IV. 



IV. SECOND PART OF THE CUTTING RULES: THERMAL FACTORS 



In the introduction wc said that both cut propagators and tic-ed propagators are on sheU, and that the difference be- 
tween them is that they carry different thermal factors. The key to determining these thermal factors is understanding 
the physical role of the cut and tic-ed propagators. We will discuss cut and tic-ed propagators separately. 



A. Cut Propagators 



Return for a moment to the one loop example discussed in the introduction. As discussed previously, the first term 
in Eqn.(Q) comes from the square of the forward scattering amplitude and represents the probability for the decay 
$ — *■ (/)i(/)2, minus the probability for the inverse decay. The field $ is the field whose self energy we are calculating. 
The fields (fii and (/)2 are the fields produced by the decay of the $ field (or the fields which combine to produce a 
$ field) and they appear as emitted lines on the right side of the scattering amplitude. This interpretation of the 
physical role played by the fields corresponding to the cut propagators allows us to understand the associated thermal 
factors. The thermal factor (1 -I- ni)(l -I- 712) — nin2 is the statistical weight associated with the probability for the 
process $ — » 0i02 minus the statistical weight for the inverse process. 

The generalization of this idea is straightforward. We write the thermal factor for the cut propagators as an equation 
of the form, 

n.Hj [(1 -I- ni)nj - n,(l + nj)] (2) 

where the product over i multiplies over fields whose momenta approach the cut line from the left, and the product 
over j multiplies over fields whose momenta approach the cut line from the right. It is easy to understand the structure 
of this expression. Performing the loop integrals puts each cut propagator on either the positive or negative mass 
shell. A particle that approaches the cut from the left and is on the positive mass shell carries a factor (1 -I- n{uj)) 
in the first term of Eqn.(|), and corresponds to an emitted particle. If the same particle were on the negative mass 
shell, the corresponding thermal factor would be (1 -t- n{—uj)) = —n{uj), which correctly reflects the fact that the 
emission of a positive energy particle is equivalent to the absorption of a negative energy paarticle. Note also that the 
momentum of a cut propagator could flow in either direction (right or left), depending on the definitions chosen for 
the momentum variables in the loops. Of course, physics must not depend on the choice of an integration variable. 
The identity 1 + n(—x) — ~n{x) ensures that all definitions of loop monenta are equivalent. 



B. Tic-ed Propagators 



Now we consider tic-ed propagators. Both cut and tic-ed propagators are on shell, and both correspond to fields 
that are emitted and absorbed. However, the fields associated with tic-ed propagators are emitted and absorbed on 
the same side of the cut self energy diagram. They represent interactions with fields in the heat bath. Naively, it 
appears that a tic-ed propagator should carry a factor N :— (1 + n) + n which would give the appropriate statistical 
weight for a field that is emitted and then subsequently absorbed. (Of course the same factor appears if emission and 
absorption occur in the opposite order). For convenience, we include a factor of 1/2 (which will be explained below) 
and use a thermal factor of the form 

N/2 (3) 

In most cases this naive guess gives the correct answer. Note that the interaction represented by the thermal factor 
N/2 is essentially trivial: the particle is produced and disappears without undergoing any interactions with the rest 
of the system. As discussed in section III, these particles are sometimes called 'spectators.' 
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There are some cases in which the fields associated with the tic-ed propagators undergo non-trivial interactions 
with the system, and in these cases the naive thermal factor is incorrect. These diagrams always have more than 
one tic-ed propagator, and thus only occur for self energies with three or more loops. The authors of H, |3) have not 
gone beyond the two loop level and thus have not seen this effect. In order to identify the cut self energy diagrams 
in which these non-trivial interactions occur, and to determine the appropriate thermal factor, we need to look at 
the corresponding scattering amplitudes. We describe below the rule for determining if the fields represented by the 
tic-ed fields undergo non-trivial interactions, and for obtaining the thermal factors for these fields. We illustrate these 
rules by looking at several cuts of the diagrams shown in Fig. [22]. 

P+Q+M 




FIG. 22: Diagrams used to illustrate the rule in section IV-B2. 



1. Notation 

We define two kinds of tic-ed propagators. 

For type-A tic-ed propagators, in the corresponding scattering amplitude, the emitted-absorbed pair of fields from 
the tic-ed propagator are 'unsplit' by either of the emitted/absorbed fields associated with the cut line. 

For type-B tic-ed propagators, in the corresponding scattering amplitude, the emitted/absorbed pair of fields from 
the tic-ed propagator are 'split' by one of the emitted/absorbed fields associated with the cut line. 

The meaning of the expressions 'split' and 'unsplit' is explained below. The cut propagators are labeled a and b 
in this section. The emitted/absorbed lines in the scattering amplitude that are associated with the cut propagators 
are horizontal arrows labeled a and b. The tic-ed propagators are labeled {x, y, z ■ ■ ■}. The corresponding emit- 
ted/absorbed fields in the scattering amplitude appear as pairs labeled {x, y, z ■ ■ ■} (see Fig. [21]). Consider a tic-ed 
propagator labeled x. In order to determine if this is a type-A or type-B tic-ed propagator, one proceeds as follows. 
Trace a path through the scattering amplitude from a to b. If you encounter both or neither of the fields labeled x, 
than the tic-ed propagator labeled x in the corresponding self energy is 'unsplit' (type-A). If you encounter one x 
field, than the tic-ed propapator is 'split' (type-B). Numerous examples of this notation are given below. 

2. The Rule 

a) If the scattering amplitude corresponds to a cut self energy that contains an arbitrary number of type-A tic-ed 
propagators, all of these propagators carry the naive thermal factor N. Several examples are shown in Fig. [23]. The 
thermal factors for these four diagrams are respectively, 

i((l + np+q)np - (1 + np)np+q)N{m)N{l -|- m) 
i((l -I- np+g)np - (1 + np)np+g)N{l)N{p + q + rn) 
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-((1 + np+q+m+k)np+m+k - (1 + np+m+k)np+q+m+k) N (p) N {k) 

^ ((1 + np+q)np - (1 + np)np+g)N{p + q + m + l + s)N{p + q + l + m)N{l) 




FIG. 23: Self energies with type-A tic-ed propagators and the corresponding scattering amplitudes 

b) I£ there is one type-B tic-ed propagator and an arbitrary number of type-A tic-ed propagators, all propagators 
carry the naive thermal factor. Several examples are shown in Fig. [24]. In the first two diagrams, the tic-ed propagator 
marked y is type-B, in the third diagram the tic-ed propagator marked x is type-B. The thermal factors for the three 
diagrams are respectively, 

(1 + np^rn-\-k))N{p)N{p + m) 
(1 + np+m+k))N{p)N{p + q + m) 
+ np+q)np - np+q{l + np))N{m)N{p + q + l)N{p + q + m + l + s) 
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FIG. 24: Self energies with one type-B ticed propagator and the corresponding scattering amplitudes. 

c) If the scattering amplitude corresponds to a cut self energy with more than one type-B tic-ed propagator, the 
thermal factors for those propagators will not have the naive form. The correct thermal factor can be written, 

2n,nj[(l + n,)nj + n,{l + n^)] (4) 

The notation is described below. Look at the set of fields coming from tic-ed propagators that appear between a and 
b. Identify the vertex in the self energy that these fields came from. The product i runs over fields whose momentum 
flows away from this vertex, and the product j runs over fields whose momentum flows towards this vertex. (Note that 
since the factor above is symmetric under interchange of i and j this definition could be reversed). This factor looks 
almost exactly like the factor for cut propagators, expect for the relative plus sign. This point and the factor of 2 in (Q) 
will be explained below. Several examples are shown in Fig. [25]. For the first diagram, the tic-ed propagators marked 
X and y are both type-B. The x in the centre of the scattering amplitude came from a vertex that the momentum 
of the tic-ed propagator flowed towards, and the y in the centre of the scattering amplitude came from a vertex that 
the momentum of the tic-ed propagator flowed away from. Thus, the thermal factor for these two propagators is 
2((1 -|- nm)''T-i+m+p+q + f^m(l + "•(+m+p+g)), SLTid the fuU thermal factor is 

((1 + np+q)np - np+q{l + np))N{p + q + m + I + s){{l + n,n)ni+„i+p+q + -|- ni+„i+p+q)) 

In the second diagram, all three tic-ed propagators are type-B. In the centre of the scattering amplitude, the prop- 
agators marked y and z have momenta that flow towards the vertex, and the propagator marked x has momentum 
that flows away from the vertex. The full thermal factor is 

2((1 + np+q)np - np+q{l + np))((l + np+q+rn+i+s)nmns + np+,+„i+;+s(l + n„i)(l + Us)) 

The thermal factors for the last four diagrams in Fig. [25] are, in order, 

2((1 + np+q+m+k)np+m+k - '^p+g+m+fe(l + np+m+fe))((l + flm)(l + "p) + nmUp) 

2((1 + '^p-\-q-\-rn-\-k )^^p+m+fc ^p+q+m+A; 

(1 + np+„+fe))((l + n,„)(l + Up+q) + Urnrip+q) 

2((1 + np+q)np - iip+q{l + np)){{l + ni)np+q+„i+i + ni{l + np+q+„i+i)) 
2((1 -t- np+q)np - np+,(l -I- np)){{l + ni)nm+i -f n;(l + Jim+i)) 
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FIG. 25: Self energies with more than one type-B tic-ed propagator and the corresponding scattering amplitudes 



3. Interpretation of the Rule 

The physical interpretation of these results is straightforward. Consider the case of two tic-ed propagators. If both 
propagators were type-A, the thermal factor would be determined by Eqn.(|^) which gives (up to numerical factors), 

N1N2 = (ei -|-ai)(e2 -faa) 

where — (1 + n^) is the statistical weight for the emission and Ui = rii is the statistical weight for absorption. 
This thermal factor tells us that type-A fields are spectators whose emission and absorption is uncorrelated with the 
emission and absorption of any of the other fields. 

If both propagators were type-B, the thermal factor would be determined by Eqn.(^ which gives (up to numerical 
factors and depending on which way the momentum is routed). 



6162 + aia2 
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This result tell us that the production of the fields associated with two type-B tic-ed propagators is correlated, and 
that the two emitted-absorbed pairs of fields interact with the system in a non trivial way. This result is consistent 
with the information we obtain from the scattering amplitude: for each of the two emitted-absorbed pairs, one of 
the partners appears on the opposite side of the emitted/absorbed line produced by a cut propagator, and thus they 
cannot be treated as spectator fields. The plus sign between the two terms corresponds to the fact that emission and 
absorption takes place on the same side of the cut line: we do not take the difference of the weighting factors for a 
given process and its inverse, but the sum of the two factors. Note that the plus sign appears in exactly the same 
way in Eqn.(|): N/2 = 1/2((1 + n) + n). 

In the case of one type-B propagator, and an arbitrary number of type-A propagators, the lone type-B field does 
not have any other field with which to correlate, and thus the factor given by Eqn.(p) is the correct one. 



Numerical Factors 



All diagrams carry an overall factor 



(27r)^+i^^ ( / dp^pidn, ] (5) 



L 



where L is the number of loops. There is an additional factor, before thermal factors are considered, which is given 

by [| 

^^|' 2^-\{^){-^Y (6) 

where p is the number of propagators and v is the number of vertices. The numerical factor contributed by the 
thermal factor is 

2{C-l)2{CorrT-l)2^ll2yjnCorrT 

where C is the number of cut propagators, CorrT is the number of correlated tic-ed propagators or the total number 
of factors involved in the products over i and j in Eqn.lQ), and UnCorrT is the number of uncorrelated tic-ed 
propagators or the number of thermal factors of the form N/2 from Eqn.(|^). The second two factors in this expression: 
2(1/2)'^"'-^°''''-^ are explicitly included in the rules given by Eqns.(||), (|^). The first two factors appear when one rewrites 
the thermal factors that one obtains directly from the Keldysh RTF Feynman rules in the physically motivated forms 
given by Eqns.(|^), (^. Combining these results and using 

C + CorrT + UnCorrT = L + 1 
L = p + 1 — V 

we find that all factors of 2 cancel and, in addition to Eqn.(|5|), we are left with a numerical factor of 

V. EXAMPLES 

In this section we give a list of the diagrams we calculated in order to deduce the rules described above. The 
calculations were performed using the Mathematica program in [^j. 

A. Two Loop Self Energy Formed by a Propagator Correction 

The two loop self energy that is formed from the one loop self energy by adding a correction to one of the propagators 
is shown in Fig. [3]. There are five cuts, as shown in Fig.[4a-e]. The cut in Fig.[4e] produces three cut propagators 
and carries a thermal factor (1 + np^g^k)nknp — rtp+g+fe(l + rip){l + n^.) according to rule in Eqn.(^. The graphs 
in Figs.[4a-d] contain two cut propagators and one tic-ed propagator. Following the rules given in Eqns. (|^), (^) the 
thermal factor for Figs. [4a, c] is Nk{{l + np+q)np — nj,+^(l -I- Up)). In the same way, the thermal factor for Figs. [4b, d] 
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is Nk+q+p{{l + np+q)np — np+g(l + Up)). For a scalar theory, the diagrams in Figs. [4a, d] can be written in the same 
form by shifting momentum integrals, and/or taking the complex conjugate. For any theory with a cubic interaction, 
these diagrams have the same structure. In the future, in order to compactify the notation, we will represent a set 
of diagrams of this form by drawing one diagram and writing a factor ^iV next to it, to indicate that there are N 
permutations of the diagram. 




(b) ( c) < J) 




FIG. 26; A two loop self energy and its cuts 



B. Two Loop Self Energy with a Vertex Correction 



The two loop self energy that is formed from the one loop self energy by adding a correction to one of the vertices is 
shown in Fig. 26a]. There are eight cuts, as shown in Fig.[26b-i]. The thermal factors are obtained from the rules given 
in Eqns.(^, Figs.[26h,i] have three cut propagators and carry the thermal factors (1 + ng+p)np-knk — np+q{l + 
rip_fc)(l + rife) and (1 + rig+fc)(l + np^k)np — Uk+qUp^ki^ + "-p), respectively. Figs.[26b-g] have two cut propagators 
and one tic-ed propagator and carry the factors: 

[(1 + nk+q)nk - nk+q{l + nk)]Np+q/2 ; 
[(1 + rik+q)nk - "-fe+g(l + nk)\Np_u/2 ; 
[(1 + nk+q)nk - nk+q{l + nk)]Np/2 ; 
[(1 + np+q)np - iip+q{l + np)]Nk/2 ; 
[(1 + np+q)np - np+q{l + np)]Nk+q/2 ; 
[(1 + np+q)np - np+q{l + np)]Np_k/'2. ■ 

In the future, we will compactify the notation by noting that Figs. [26b, d,e,f], and Figs. [26c, g] are permutations of 
each other. We will draw Fig. [26b], with a factor ®4 and Fig. [26c], with a factor Cg)2. 



C. Further Examples 



The thermal factors for the diagrams listed below are determined from the rules given by Eqns.(||), (||) and (^. To 
simplify notation, we define the symbol 

J\f{xi,- ■ -XnWi, ■ ■■Vm) = [(1 + n{xi) • • • (1 + x{x.n)n{yi) ■ ■ - yixn) - n{xi) ■ ■ •n(a;„)(l + n{yi)) • • • (1 + n{y„i))] 



'(b) ' (d) 

FIG. 27: A three loop self energy and its cuts. 

The thermal factors for Figs.[27b-d] are respectively, 

-N{k)M{p + q + m;m,p) 
^N{k)N{m)M{p + q;p) 
-^N{m)N{k)N{p + q;p) 




FIG. 28: A three loop self energy and its cuts. 

The thermal factors for Figs.[28b-e] are respectively, 

M{p + q + m,l;p,l + m) 

-A^(p + q + m)M{p + q + I + m\p,l + m) 

^N{m)N{l + m)M{p + q;p) 

^N{l)N{m)U{p + q;p) 
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P+M+Q+L 




(c ) 



FIG. 29: A three loop self energy and its cuts. 



The thermal factors for Figs.[29b-d] are respectively, 



M{p + q + m + l;p,m,l) 
^N{l)N{p)N'{p + q + m;p + ni) 



N{l)JV{p + q + m;p,m) 




P+Q 




(b) 



(e) 




( c) 






(d) 




(f) 



(g) 



FIG. 30: A three loop self energy and its cuts. 



The thermal factors for Figs. [30b- g] are respectively, 

A/'(p + q,m,k;p + m + k) 
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^iV(fc)A/'(p + q,m;p + m) 

jN{k)N{m)N'{p + q;p) 

^N{m)N{p + m + k)Af{p + q;p) 

^N{k)N{p + m)M{p + q;p) 

2[(1 + nk)np+rn+k +nk{l + np+rn+k)W{p + q;P) 




FIG. 31: A three loop self energy and its cuts. 

The thermal factors for Figs. [31b- i] are respectively, 

J^{p + q + m + k;k,m,p) 
^N{p)M{p + q + m + k;k,p + m) 

^N{q + p)N{q + p + m + k)M{p + q + m;p + m) 

^N{m)N{k)N'{p + q + m + k;p + m + k) 

^N{p + q)N{p + q + m)N{p + q + m + k;p + m + k) 

^N{p)N{k)N'{p + q + m + k;p + m + k) 

^N{m)N{p + q + m)N{p + q + m + k;p + m + k) 

2[(1 + nm){^ + rip) + nmnp]N{p + q + m + k;p + m + k) 







FIG. 32: A four loop self energy and its cuts. 
The thermal factors for Figs.[32b-ni] are respectively, 

^N{k)N{p + m)N{v + q;p) 
^N{k)Af{p + q,m;p + m) 

^N{p + m + q)M{p + q + k;p,k) 
J^{p + q + k,m;p + m,k) 
^N{k)N{m)N'{p + q;p) 

^N{k)J\f{p + q,m;p + m) 

^N{k)N'{p + q + m;p,m) 

^N{k)N{m)J\f{p + q + m;p + m) 
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^N{k)N{p)Af{p + q + m;p + m) 

jN{k)N{p + q)M{p + q + m;p + m) 

jN{k)N{p + q)M{p + q + m;p + m) 

2[(1 + nfc)n,+p+fe + nfc(l + n,+j,+fe)]A/'(p + q + m;p + m) 




FIG. 33: A four loop self energy and its cuts. 

The thermal factors for Figs.[33b-p] arc respectively, 
J\f{p + q + m + I + s; s,m,l,p) 
^N{s)N'{p + q + m + l;m,l,p) 

2[(1 + ns)np+g+m+i+s + ns{l + np+q+m+i+s)W{p + q + l;l,p) 
^N{p + q + m + l + s)N{m)Af{p + q + l;l,p) 

^N{s)N{m)Af{p + q + l;l,p) 

^N{p + q + m + l + s)N{p + q + l + m)N{p + q + I] l,p) 

{l+ns){l + nm)W{p + q;p) 
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In{p + q + m + l + s)N{m)N{l)J\f{p + q;p) 
8 

^Nis)Nim)Nil)N-ip + q;p) 

In{p + q + m + l + s)N{p + q + m + l)N{p + q + l)N{p + q]p) 
o 

\n{p + q + m + l + s)N{p + q + m + l)N{l)N{p + q\ p) 
o 

(1 + nrr,)\N{p + q;p) 

lN{m)N{p + q + m + l + s)N{p + q + l)M{p + q\ p) 
8 

N{p + q + + np+q+„i+i+s)ns + ?ip+g+„i+;+s(l + ns)]Af{p + q;p) 



VI. CONCLUSIONS 



In this paper we have discussed a set of rules for calculating the imaginary parts of self energy diagrams as a series 
of tree amplitudes which represent physical scattering processes. The thermal factors associated with each scattering 
diagram have a physical interpretation in terms of the statistical weighting factors associated with the emission and 
absorption of thermal fields. Work on a general derivation of these rules from first principles is in progress. 
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